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Canonical ensemble of an interacting Bose gas: 
stochastic matter fields and their coherence 
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We present a novel quantum stochastic evolution equation for a matter field describing the canoni- 
cal state of a weakly interacting ultracold Bose gas. In the ideal gas limit our approach is exact. This 
numerically very stable equation suppresses high-energy fluctuations exponentially, which enables 
us to describe condensed and thermal atoms within the same formalism. We present applications 
to ground state occupation and fluctuations, density profile of ground state and thermal cloud, and 
ground state number statistics. Our main aim are spatial coherence properties which we investigate 
through the determination of interference contrast and spatial density correlations. Parameters are 
taken from actual experiments Q]. 
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Equilibrium fluctuations in ultracold gases reveal de- 
tailed information about states and phases of interacting 
many-body quantum systems Q. Recent experiments 
permit to control ultracold quantum gases in a hitherto 
unknown precision and to investigate temperature depen- 
dent quantities like the thermal density and ground state 
occupancy [3L spatial correlation functions (3-0) density 
fluctuations TQ or interference contrast 

In this work we determine the canonical state of an in- 
teracting ultracold Bose gas. A novel quantum stochas- 
tic evolution equation for a c- number field ip(x) is pre- 
sented such that canonical quantum statistical expecta- 
tion values can be replaced by an ensemble mean over 
these stochastic matter fields. The equation allows to 
determine coherence properties and other relevant ob- 
servables; it is based on a mean-field type approximation 
and strictly valid for the non-interacting case. 

Most theoretical descriptions of interacting ultracold 
Bose gases at finite temperature are based on grand 
canonical statistics @, [t| . For actual experiments involv- 
ing a fixed and finite number of particles, however, a 
canonical description is natural. In studying the role of 
the chosen ensemble, attention so far has been paid to 
ground state number fluctuations [l(| E[ ■ While for the 
ideal gas canonical and grand canonical ensemble give 
vastly different predictions [IH , this ceases to be true for 
interacting gases in the thermodynamical limit [i~3| . Our 
work is based on canonical statistics right from the start 
and allows us to not only investigate occupation fluctua- 
tions but also spatial coherence properties. 

Many stochastic field methods exist for the description 
of Bose gases at finite temperature. All of these are based 
on grand canonical statistics, and nicely overviewed and 
compared in [j| [l4[ . In the truncated Wigner approach, 
the evolution of the field Wigner functional is determined 
approximately from a sampling over random initial fields 
whose dynamics is given by the Gross-Pitaevskii equation 
[l5l [l6j . Based on a quantum kinetic theory and a separa- 
tion of condensed and non-condensed part, Gardiner and 
co-workers derive a stochastic Gross-Pitaevskii equation 
fl ij ■ With a similar result, a functional integral approach 



to the evolution of the field Wigner distribution js worked 
out in an approach by Stoof and co-workers (Tsj . Care 
has to be taken with respect to the white noise driving 
these equations. Exact methods based on the positive 
P-representation are used by Drummond and co-workers 
19]. It is possible to use this approach for 3D systems; 
still, the long-time numerical solution has to be exercised 
with caution. We see the strength of our approach in its 
unified applicability to a vast number of different phe- 
nomena: from ground state fluctuations to properties of 
the thermal cloud, to coherence properties and contrast 
in Bose gas interferometry. For the latter we obtain nice 
agreement with experiments of the Schmiedmayer group 
that may be well described by Luttinger liquid theory [l| 
or by a stochastic phase model [20]. 

Two properties of our novel equation should be empha- 
sized: first, unlike in our previous attempt 21], the equa- 
tion is not norm preserving. Still, the norm fluctuations 
are small compared to those of related stochastic equa- 
tions used for grand canonical simulations. Secondly, as 
in [21], ultraviolet cutoff problems do not appear due to 
the use of the Glauber-Sudarshan P-function: effectively, 
our treatment leads to spatially correlated noise, unphys- 
ically large momentum kicks are suppressed. These prop- 
erties afford a very stable numerical solution of full 3D 
problems, using arbitrary trap potentials. Due to lack 
of space and its current interest, we here concentrate on 
ID gas interference. Still, we emphasize that we are also 
able to treat full 3D gases within our approach [22| . 

We propose to use the stochastic (Ito) equation 
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for a c-number matter field ij)(x,£) = (x\ip(t)} to deter- 
mine all equilibrium properties of a weakly interacting 
Bose gas of N particles at arbitrary temperature T (k is 
Boltzmann's constant). Throughout, we will refer to ([1]) 
as the stochastic matter field equation (SMFE) for finite 
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temperature. Crucially, the operator H is the effective 
(mean-field) one-particle energy operator 



H=7^ + V(x)+g(N-l) 



(m\m) 
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such that equ. (JlJ may also be seen as a stochastic Gross- 
Pitaevskii equation. As usual, V{x) denotes the trap po- 
tential, the interaction parameter g is proportional to the 
s-wave scattering length a s and m is the mass of a Boson. 
The parameter A appearing in (fT]) is a phenomenological 
damping rate that sets the time scale for transition to 
equilibrium. Its appearance as square root with the fluc- 
tuations reflects a fluctuation-dissipation-relation. The 
fluctuating part is driven by complex Ito increments \d^) 
with \d£){d\\ = Idt, | = 0. Note, however that 

the operator %/ He~ H / kT acts on the noise, effectively 
leading to spatially correlated noise (2lj . 

Before we show the versatility and accuracy of the 
SMFE in applications later, let us sketch how we ar- 
rive at (fTJ). Our aim is to determine mean values 
(. . .) jv =tr[. . . pn] with the canonical density operator 
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in second quantization with Hamiltonian H, canon- 
ical partition function Zjy, and projector II jv = 
J2 l{ n fe})({ n fe}l on t° the iV-particle subspace. 

Normally-ordered matter field correlation functions are 
expressed in terms of functional phase space integrals 
21], for instance 
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and Cn = J T^W\ Wn{4>)- Note that second (or higher) 
order correlations require the use of Wn-2 (or lower 
index) in expression while Cm remains. 

For the ideal gas case (g = 0), we prove that the SMFE 
([T|) corresponds to a Fokker-Planck equation [24j whose 
stationary solution is just the weight functional Wjv(V')- 
Thus, equilibrium expectation values of the canonical en- 
semble are obtained from propagating equ. ([1]) and av- 
eraging. In practice, we use a long-time-average over a 
single trajectory ip(x,t). 

The SMFE (JTJ) is exact for an ideal gas; interactions 
can be included with great success: we propose to use 
the single equation (TTJ) with ©, containing the current 
stochastic mean field energy, to describe all properties 
of weakly interacting Bose gases in a unified way. In- 
deed, we emphasize that (TTJ interpolates smoothly be- 
tween the high-temperature limit T » T c , when inter- 
actions are negligible and thus our description is exact 
anyway. At the opposite end, when T <C T c , the SMFE 



(fTJ) with ([2]) is just mean field Gross-Pitaevskii theory 
and is again expected to give good results. Clearly, 
the fluctuations we describe are of thermal origin; quan- 
tum fluctuations are taken into account to some extent 
through the use of the P-representation. Note also that 
it is of crucial importance to use the current, stochas- 
tic if>(x,t) in @, such that on average (i?)(jv-i) = 
p 2 /2m+V(x)+g(ilj^(x)'ilj(x)}( N _iy The quality of choice 
@ was tested by solving equation (TTJ) for a two mode 
system and comparing with numerically exact quantum 
results over a wide temperature range. 
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FIG. 1: Ground state occupation (top, left), its variance 
(top, right) Ano = y ((no — (no)) 2 }, third (bottom, left) and 
fourth (bottom, right) central moments as a function of tem- 
perature (scaled with the critical temperature of the thermo- 
dynamic limit T c = hu)N 1/3 /k((3) 1/3 ) for an ideal 3D Bose 
gas of 200 particles in a harmonic trap. The data obtained 
with the SMFE (black plus signs) is compared with the results 
of the recursion relation (red solid line). 

We convince ourselves of the validity of ([T]) by first 
considering an ideal Bose gas of 200 particles in a 3D 
harmonic trap. In Fig. Q] results for ground state oc- 
cupation, its variance and further centered moments are 
compared with exact results for the canonical ensemble 
obtained from a recursion relation (25j . 

As a first application to the interacting case in Fig. 
[2 the density profile (green solid line) of a 87 Rb quasi- 
1D gas of 20240 atoms in a trap with frequencies oj z — 
2tt x 9Hz and ui± = 2ir x 32Hz at a temperature of 
185nA is shown (we use ID coupling constant gm = 
2hco±a s ). The blue dashed line is the contribution with 
off-diagonal long range order (ODLRO) whose wave func- 
tion ipo( z ) is obtained from a diagonalization of the full 
p(z,z') = (ibHz)il){z'))N applying the Penrose-Onsager 
criterion [26||. Moreover, a histogram of the stochastic 
occupation n = (ipo\p\ipo) leads to the ground state 
number statistics P(jiq) (inset in Fig. [2|. The average 
number of particles in state ipo turns out to be 12573. 
Our findings are nicely com pat ible with the "stochastic 
Gross-Pitaevskii" results of [14j |. without, however over- 
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FIG. 2: Density profile (green solid line) of an ID interacting 
Bose gas of 87 Rb atoms with uj z = 2-7T x 9Hz,tv± — 2-k x 3QHz 
at a temperature of 185nK simulated with our SMFE. The 
ground state contribution (blue dashed line) and the thermal 
part (red dashed-dotted line) are obtained with the Penrose- 
Onsager criterion [26j] by calculating the full density matrix 
p(z,z'). In the upper left corner the ground state number 
statistics P(no) from our simulation is shown. 
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FIG. 3: Length dependence of the average contrast (|j4(L)| 2 ) 
of an interference pattern of two uncoupled ID Bose gases in 
a harmonic trap for temperatures of 33nK, 47nA" and 68nK 
(niD ~ 50 fim~ 1 , u}± — 2k x 3.0kHz). The data from the 
SMFE (black plus signs, green crosses, brown stars) is com- 
pared to Luttinger-liquid theory (LLT) (red solid line, blue 
dashed line, yellow dashed-dotted line) which agrees well with 
the experimental measurements [lj. Small deviations arise 
from the variation of the density in the harmonic trap, while 
the Luttinger-liquid theory applies to a uniform density. 



estimating lowly occupied regions. 

The SMFE ([1]) is ideally suited to study coherence 
properties of interacting matter waves through the de- 
termination of spatial correlation functions. As an ap- 
plication we show results of our SMFE applied to re- 
cent experiments in the Schmiedmayer group [1]: two 



FIG. 4: Distribution functions of the interference contrast for 
different lengths L and different temperatures. The length- 



dependent normalized interference contrast a 



(\W) 



■zr is sam- 



pled over 10000 realizations of the SMFE. The calculation 
is done for different temperatures and different integration 
lengths L (nm = 59/wn -1 , uj± — 2tv x 3.0kHz). 



independent condensates are prepared in quasi- ID; af- 
ter expansion they interfere; the observed interference 
pattern is integrated over a length L which determines 

L/2 

the contrast |A(L)| 2 , where A(L) = J dzip 1 (z)tl)2(z). 

-L/2 

Both, mean value (\A(L)\ 2 ) (Fig. [3]) and the full distribu- 
tion VF(a) of the normalized moments defined through 



/ W(a)a m da 



<l^l 2 >" 



are determined (Fig. 



2]). In [l| it is shown that experimental results are well 
described by Luttinger-liquid theory (LLT) to which we 
will compare our SMFE results. 

In Fig. [3] we show the average contrast (|t4.(L)| 2 ) as a 
function of length L for different temperatures and find 
very good agreement with LLT (and thus with experi- 
ment). Deviations for large L arise from density varia- 
tions along the gas: LLT results are based on a uniform 
density. The gas contains some 4400 87 Rb atoms with 
a central density of nm ~ 50/ito -1 . The interaction 
strength for this ID case is again giu — 2hujj_a s , with 
ui± = 2tt x 3.0kHz. In Fig. 2] we show the distribution 
function of the interference contrast VF(a) as a histogram 
with a{L) = \A(L)\ 2 / (\A(L)\ 2 ) obtained from our sim- 
ulations. We use a central density riic = 59/im -1 in 
line with the experimental setup (see [l|). LLT predicts 
a change of the shape of the distribution function with 
decreasing parameter F = which is excellently 

reproduced by our simulations. 

Finally, as in 0, we investigate two-point den- 
sity correlation functions after expansion g2(x;t) — 

T^P^^W^ - We chose * - 27ms; note that 
free (non-interacting) expansion can be assumed. A gas 
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FIG. 5: Normalized spectrum of density ripples 
i\p(<l)\ 2 } /( n iD^h) for a weakly interacting Bose gas ( 87 Rb) 
with a central density of mo w 40/im -1 , in a trap with 
transversal frequency aii = 2n x 2kHz; the expansion time 
is t = 27ms, healing length is £h = \fl^- The simulations of 
the SMFE are done for 12nK (black plus signs), 27nK(red 
crosses) and 40nK (blue stars) and compared to theory of 
[27l | (YlnK blue solid line, 27nK black dashed line and 40^ 
red dashed dotted line). 



of about 10000 atoms in a harmonic trap with wj_ = 
2ir x 2.0kHz and a central density of nir> w 40^m~ 1 is 



used. As in [23], we Fourier transform to obtain den- 

00 

sity ripples (\p(q) 2 \) := n\ D J exp(iqx)[g 2 (x;t) - 1]. In 

—00 

Fig. [5]our results are compared to calculations from [27]. 
Again, we see very good agreement. 

Let us summarize our achievement: We present a novel 
quantum stochastic matter field equation (SMFE) for a 
gas of N particles trapped in an arbitrary potential at 
any temperature T (canonical ensemble). The equation 
is strictly valid for a non-interacting gas; we include in- 
teractions in a stochastic mean-field sense and obtain 
promising results over the entire relevant temperature 
regime. The SMFE is capable of tackling problems in 
ID to 3D with arbitrary trapping potentials. Results for 
ground state occupation distribution and density profiles 
are shown. Of particular interest is the determination of 
spatial correlation functions of arbitrary order. We apply 
our approach to calculate interference contrast distribu- 
tions and density ripples as recently measured. Our re- 
sults are in good agreement with Luttinger liquid theory 
which has proven to describe the experiments adequately 
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